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abstract 

We study the effects of an integrable impurity in a periodic t-J cliain. The impurity 
couples to both spin and charge degrees of freedom and has the interesting feature that the 
interaction with the bulk can be varied continuously without losing integrability. We first 
consider ground state properties close to half-filling in the presence of a small bulk magnetic 
field. We calculate the impurity contributions to the (zero temperature) susceptibilities 
and the low temperature specific heat and determine the high-temperature characteristics 
of the impurity. We then investigate transport properties by computing the spin and charge 
stiffnesses at zero temperature. Finally the impurity phase-shifts are calculated and the 
existence of an impurity bound state in the holon sector is established. 
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1 Introduction 



Theoretical investigations of strongly correlated electron systems have shown that the low temperature 
properties of such one dimensional systems have to be described in terms of a Luttinger liquid rather 
than a Fermi liquid. Of particular interest also from an experimental point of view are the transport 
properties of these systems in the presence of boundaries and potential scatterers. Several attempts 
have been made to describe such a situation: the transport properties of a ID interacting electron gas 
in the presence of a potential barrier have been first studied using renormalization group techniques 
[Q, ^ . Triggered by a more recent study of this problem by Kane and Fisher different approaches 
such as boundary conformal field theory and an exact solution by means of a mapping to the 
boundary sine-Gordon model ^, U, have been applied to this problem. In particular the low 
temperature properties of magnetic (Kondo) impurities in a Luttinger liquid |9|, |l^, 12 1 have been 



investigated in great detail. In the present work we will investigate the effects of a particular type of 
potential impurity in a Luttinger liquid (where both spin-and charge degrees of freedom are gapless) 
by means of an exact solution through the Quantum Inverse Scattering Method (QISM). Attempts 
to study effects due to the presence of impurities in many-body quantum system in the framework 
of integrable models have a long successful history 14, 15, [l^, 0, ITs). As far as lattice models 



are concerned the basic mechanisms underlying these exact solutions are based on the fact that the 
QISM allows for the introduction of certain "inhomogeneities" into vertex models without spoiling 
integrability. Two approaches are possible and have been studied for various models: 
Since the /^-operators defining the local vertices satisfy a Yang-Baxter equation with an TZ matrix 
depending on the difference of spectral parameters only, one can build families for vertex models with 
site-dependent shifts of the spectral parameters. This has been widely used in solving models for 
particles with an internal degree of freedom by means of the nested Bethe Ansatz [19] and, more 



recently, also for the construction of systems with integrable impurities [16| (for a particularly simple 
case see ||20t| ). 

A second possibility is to introduce an impurity by choosing an ^-operator intertwining between 
different representations of the underlying algebra (possibly with an additional shift of the spectral 
parameter). This has first been used by Andrei and Johannesson to study the effect of a spin S site in 



a spin-i Heisenberg chain |13|. Later this approach was used to study chains with alternating spins 

[0, m ^ 

Our work is based on the second approach. A novel feature as compared to |l^ is the presence of 
a continuously varying free parameter describing the strength of the coupling of the impurity to the 
host chain. The presence of the free parameter is based on properties of the underlying symmetry 
algebra of the model, which in our case is the graded Lie algebra gl{2\l) (see below). 

Recently, vertex models and the corresponding quantum chains invariant under the action of 
graded Lie algebras have attracted considerable interest 21, |22|, p^, p3, [2^, 27]. In addition to 



the well known supersymmetric t-J model [| ]|2^, 3C, 31, 32, |33| (which is based on the fundamental 
three-dimensional representation of gl {2\1)) a model for electrons with correlated hopping has been 
constructed using the one-parametric family of four-dimensional representations of g/(2]l) ]p3, |3^, Rl, 



In this paper we study the properties of the supersymmetric t-J model with one vertex replaced 
by an C operator acting on a four-dimensional quantum space. This preserves the supersymmetry of 
the model but at the same time lifts the restriction of no double occupancy present in the t-J model 
at the impurity site. The free parameter associated with the four-dimensional representation of the 
superalgebra allows one to tune the coupling of the impurity to the host chain. Note that the present 
models allows for the study a situation which in some respects is more general than the ones mentioned 
above: the impurity introduced here couples to botli spin- and charge degrees of freedom of the bulk 
Luttinger liquid. The extension of our calculation to the case of many impurities is straightforward 
(and will in fact be used in some parts of the paper). 

The paper is organized as follows: In the following section we give a brief review of the Quantum 
Inverse Scattering Method and the Bethe Ansatz for 51/(2] l)-invariant models. Then we study the 
ground state properties, low temperature specific heat and transport properties and show how they 
are affected by the presence of the impurity. Finally we compute the phase shifts acquired by the 
elementary excitations, holons and spinous, when scattered off the impurity. 



^For a collection of reprints see jESl 
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2 Construction of the model 



The impurity model is constructed by means of the graded version of the Quantum Inverse Scattering 
Method ||3^, m 41]. We start with the 7^-matrix of the super symmetric t-J model [32] 



ntj = n. 



33 



A + i 



n + 



X + i' 



(2.1) 



with n^^^^ = 5aif)2<^a26i (— l)'^''i'^''2 being the graded permutation operator acting in the tensor product of 
two "matrix spaces" isomorphic to the three-dimensional "quantum space" space spanned by (TiiiO) 
with the respective grading = = 1 and Eq = 0. The corresponding t-J ^-operator is given by 
^33 = 117^33 and satisfies the following (graded) intertwining relation 



7^33(A - fi) 





^33 (A - fi) 



In components this equations reads 

7^33(A-M)^^S^33(A)^V''£33(/^)^f (-l)^-(^-+^^^) = /:33(/^)?,'^^/:33(A) 



a.2C2 I 
7/3 ( 



C262 ■ 

(2.2) 

The monodromy matrix Ttj is defined on the graded tensor product on quantum spaces and its 
matrix elements are given by 

'TtA^r.Zl' = ^33(A)r/,,£33(A)r_-^^_^ • • • £33(A)-^,^(-l)^'-K+^^^) T^l . (2.3) 

The hamiltonian of the t-J model is then given as the logarithmic derivative of the transfer matrix 
r(A) = str[Ttj(A)] := ELiC-I)"" ['^ij(A)]'"' at zero spectral parameter Jpl 



n 



tj 



-.|^ln(r33: 



A=0 



2AA, 



J2{^k,k+l - I) - 2Me 



(2.4) 



k=l 



The 7^33-matrix can be constructed as the intertwiner of the three-dimensional fundamental repre- 
sentation of the superalgebra gl{2\l). An impurity model can then be constructed in a way analo- 
gous to the spin-S* impurity in a spin-i Heisenberg chain ]|l^ by considering the intertwiner of the 
three-dimensional representation with ithe typical four-dimensional representation corresponding to 
an impurity site with four possible states (T)i,Ti= 2,0). In this way we obtain an ^-operator £,34 
with the same auxiliary space dimension as the 7?.33-matrix satisfying the following equation: 



^3a{X) 



^33(A-/i) 




C34{X) 




^33 (A - /i) 



The double line denotes the four-dimensional space (t, |, 2, 0) with the respective grading = ei = 1 
and £2 = £0 = 0. So double occupancy is possible at the impurity. The i234-matrix is given by 



34 



A-i(f + 1) 

A + i(f + 1) ' A + i(f + 1 



+ 



(2.5) 
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where £, expressed in terms of projection operators (the so called 'Hubbard Operators') X"-^ = \a){b\ 
with o, b =t, I, 2, 0, is given by 



C 



-4^ 



.10 



V^TTx!^^ - ^xf \ 
^xf + V^TTx°-^ 



(2.6) 



\ \/aTTxf - ^^2^ + VoTIXf a + + X\^ + 2Xf / 



The parameter a > is associated with the four-dimensional representation of (7/(2|l)|^ By 
inserting an £34-operator at one site (for example site 2) we arrive at the following monodromy matrix 
for the impurity model (see Fig. ^ 

%^vi^)Z::::;:tl = ^^^^^^^ 

(2.7) 

The hamiltonian is given by the logarithmic derivative of the transfer matrix at zero spectral parameter 



n 



tj~lmp 



.ain(r(A)), 



(2.8) 



Due to the fact that £34 has no shift-point the hamiltonian contains three-site interactions of the 
impurity with the two neighbouring sites. This makes the local hamiltonian rather complicated. 
However, as is shown in Appendix C by direct computation, the hamiltonian constructed in the way 
described above is invariant under the graded Lie-algebra gl{2\l). Also in the continuum limit only a 
comparably small number of terms relevant in the renormalization group sense will survive. 
The impurity contributions to the hamiltonian are found to be 



' - (^3-i(o))::>s3 (-i)[(-3-.3)^.^] - . (£3-/(0));;;:^ {c',M)it 



(2.9) 



which can be expressed in terms of the Hubbard operators as 

a + l 1 



-Hi 



,2,3 



1 + 



(f + 1)' (f + 1) 



ni2 + aX?o + Hu + ^{l- 2X™) ni3(l - 2X?0) 

+ (n23 + x|2(i + ni3)) - f [(1 - 2X°o) ni3ni2 + h.c. 



a 
2 



1 - 2X?° j ni3i^i2 + h.c.\ + HuUrsHu + [H12U13U12 + h.c. 
Here ni2 denotes a modified permutation operator 

ni2 = + (1 - x|2)ni2(i - x|2)) . 



The operator H12 is given by 

H12 = xl°\{V^^-i)x^^ 



xi" iV^TT - i)x'^^ + V^x. 



42 



h.c. 



(2.10) 



(2.11) 



From the explicit form of the hamiltonian it is clear that the impurity can be thought of as being 
attached from the outside to the t-J chain as is depicted in Fig. |l| 

In the limiting cases a ^ 00 and a ^ the impurity contribution simplifies essentially 



lim i^i 2,3 
lim 2,3 



(1 - n2,tn2,|)(2 - ni2 - n23)(l - "2,1^2,1) + n2,tn2,j,(3 - ni3), 
l-(l-2X?°)ni3(l-2X?°). 



(2.12) 



Other parameter regions may be possible, but will not be considered here. 
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It is shown below that for q = the impurity site behaves hke an "ordinary" t-J site in the ground 
state below half-filling. In the limit a = oo the impurity becomes doubly occupied and causes the 
hopping amplitude between the two neighbouring t-J sites to switch sign. The situation is equivalent 
to a decoupled impurity and & t-J chain with L sites and a twist of —1 in the boundary conditions. 
These two situations are shown in Fig. ^ and Fig. ^ respectively. 

For the half-filled case (two electrons at the impurity site) the impurity decouples from the host 
chain and we obtain the following hamiltonian (see Fig. |^ 

^^fer"^ = -ni3 + i + ^, (2.13) 

If we wish to consider the case of many impurities the above construction will go through as long 
as every ^-operator VC34 is sandwiched by two vC-operators £33. The resulting hamiltonian will consist 
of a sum over terms of the form Hk^k+i,k+2 and 1 — Iljj+i. The maximal number of impurities is thus 
bounded by the number of t-J sites. 



3 Algebraic Bethe Ansatz 



The Bethe Ansatz equations can be derived as in [32|. The monodromy matrix is a 3 x 3 matrix of 
quantum operators acting on the entire chain 



imp 



( ^11 (A) 



A22{X) B2{X) 

C2{X) D{X) I 



(3.1) 



Like in the case of the t-J model one constructs the eigenstates of the hamiltonian by successive 
application of the Cj-operators on the bosonic reference state |0) 



|Ai, . . . , A„|F) = CaAXi)CaA><2) ■ ■ ■ C7a„(A„)|0)F'^- 



■ ai 



(3.2) 



The values of the spectral parameters Aj and the coefficients j^^-^-.-ai^ which are constructed by means 
of a nested algebraic Bethe ansatz (NAB A), are determined by requiring the cancellation of the so- 
called "unwanted terms" . For the construction of the NABA one only needs the intertwining relation 
( |2.2D (and the existence of the reference state). Due to the fact that the corresponding 7^-matrix is 
the same for the impurity model and the t-J model the resulting equations are very similar. The 
only modification of the Bethe equations arises from a different eigenvalue of the ^jj-operators on the 
reference state leading to the following form for the Bethe equations 




n 



A. 



Aa 



Ne x(l) 

n ° 
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=1 Xj 

N 

■ n 



A« + 



i,...,Ne = m + N^ 



Aa 



Xf+i 



a 



The energy of a Bethe state of the form (3.2) with spectral parameters {A 
1 . . . A'^l } in the grand canonical ensemble is given by 



'31 



1 '^i + i 



iVi- (3.3) 
i = l...N,}, {A^') ,7 = 

(3.4) 



Here /x is the chemical potential and H is a bulk magnetic field. Properties of the ground state and 
excitations as well as the thermodymanics can now as usual be determined via an analysis of the Bethe 
equations (^^). In order to analyze ground state and excitations we first have to discuss some details 
concerning the lattice length. We note that the length of the lattice is L + 1, where L is the length of 
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the t-J "host" chain. It is known that the ground state of the t-J model changes if the length of the 
lattice is increased by one. The situation is similar in the present model. A unique antiferromagnetic 
ground state exists for odd (equal) numbers of up and down spins. If we require a smooth limit to the 
half-filled band with a doubly occupied impurity we find that the length of the host chain must be of 
the form L = mod 4. 

If we consider the case of Ui impurities and t-J sites discussed above the left-hand-side of the 

first set of Bethe equations changes to — f- I — I whereas the other equations remain 

\-^o + 2 J V — 2~V 

unchanged. 



4 Ground State Properties 

In order to study ground state properties we need to know the configuration of Aj's and A^^^'s corre- 
sponding to the lowest energy state for given /i and H. It can be found in complete analogy to the t-J 
model without impurity: for finite magnetic field the ground state is described in terms of two filled 
Fermi seas of 

• complex A-A(^)- "strings" A± = A^^^ii fsO] (where A^^) are real and fill all vacancies between the 
Fermi rapidity B and oo and — oo and —B). As we approach half- filling B tends to 0. Removing 
one A-A^^)- "string" from the Fermi sea and placing it outside leads to a particle-hole excitation 
involving only charge degrees of freedom ( "holon-antiholon" excitation)[31]. 

• real solutions Xj associated with the spin degrees of freedom. They are filling all vacancies 
between A and oo and — oo and —A. We note that A — > oo as — > 0. 

At zero temperature the dressed energies ^'(A) and e(A) of the excitations associated with charge and 
spin degrees of freedom respectively have to satisfy the following coupled integral equations pO| ] 

rao p — B roc r — A 

'I'(A) = 27ra2(A) - 2/i - / +/ a2{X - n)'^{ti) - + ai(A - ^)e(/u), 

H fOor-B 

e(A) = 27rai(A)-^--- /+/ ai(A - /z)^'(/u) (4.1) 

2 Jb J-oo 

where a„(A) = ^ "'^■^ . The functions ^' and e are negative in the intervals {—oo,B) U {B,oo) and 

(— oo,^) U (A, oo) as they are monotonically decreasing functions of |A|. 
The ground state energy per site is given by 

E - fiNe- HS"" , if 2a 



L 



If 2a \ 

^'(0)-2^ + 2 + - -2;U + 2 — -/ d^a„(^)^'(^) . (4.2) 

L \ a + 2 J-B 



The first part on the r.h.s. is the contribution of the L sites of the host t-J chain, whereas the 
second term is the contribution from the impurity. Transforming the integral equations into the 
complementary integration intervals we obtain the following equations 



/A rB 
dfjL a2{\ - /u) Esili) + / dfjL ai(A - /i) ec(/i) , 
^A J-B 

H 

ec{\) = — + d^i ai{X- fi) Esin) , (4.3) 
2 J-A 

where Eg = —e and Ec = The ground state energy per site is now given by 

= e.(0)-2, + 2 + i(.„ + -l^-2,) . (4.4) 
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where 



B 



-B 



(4.5) 



We now split ( [4. 4]) into two parts: the contribution of the t-J host chain Chost and the contribution of 
the impurity Eimp. Note that in ( [4.4D we divided by the length of the host chain instead of the length 
of the lattice. The impurity contribution to the ground state energy is then given by 



imp 



a + 2 



2/i 



(4.6) 



Before we turn to an analysis of the impurity contributions to the ground state energy we give a 
brief review of the properties of the t-J host chain. The ground state properties of the t-J model have 
been studied in great detail in ||4^, 44]: below a critical density nc, which is related to the magnetic 
field by 



H = 4sin^ 



irrir 



(4.7) 



the integration boundary B in (^^) is oo and the system behaves like a Fermi-liquid with all spins 
up. The particle density is a function of A in this case. For densities larger than ric the system is a 
Luttinger liquid exhibiting spin and charge separation. For general band-fillings the integral equations 
( ^.31 ) can be solved only numerically (in the almost empty band it is possible to reduce them to a 
system of coupled Wiener-Hopf equations but we do not pursue this avenue here) . For densities slightly 
abov e ric {B ^ 1) and B ^ A the integral equations can be solved analytically. Evaluating equation 
(O) in this case leads to: 



ec(A) 



-27rai(A) + H + i 



H 

y 



1 



1 

Vbj 



+ o 



1 

l3 



H 



d^i ai(A - /i) Esili) , 



(4.8) 



Via Ssi^) = and Sc{B) = the integration boundaries depend on the magnetic field and the chemical 
potential in the following way 



Ar 



^/4^rg 



+ ••• 



Br 



4arctan(2ylc) - \/^\/4 - H 



2TTy^fl - H/2 



+ 



(4.9) 



Using the densities pc and ps (see ( 4.23 )) the zero-temperature density and magnetization per site are 
given by: 



D, 



bulk 



1 



1 



1 



irBr 



vr 



■ arctan(2Ac)) 



mbuik 



1 + 



1 \ 1 



TlBr 



vr 



■ arctan(2Ac 



(4.10) 



oo m 



Taking the limit H ^ with finite B^ is not permitted in ( [4.10D as this would imply A 
contradiction to B ^ J^. 

For the almost half-filled band and a small magnetic field ( |4.3| ) can be solved analytically as well 
[^ ]: combining a Wiener-Hopf analysis for small magnetic field (^ ^ 1) and an iterative solution 
near half-filling (B <^ 1) one obtains |4(:] 



where p, = 21n(2) — p, \p\ <^ 1 and 



e 



(1 



1 



2ln{H) 



+ ...) 



-p — 2a 



B' 



21n(2)^53 



3C(3) 



p + 2a + 



81n(2) 



3^ V6C(3) 



{p + 2ay- 



(4.11) 



(4.12) 



'The correct result would be D^uik ~ and not Dtuik = (see 
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From these results we can now determine bulk and impurity contributions to the ground state energy. 
In order to give a reference frame for the impurity properties we start by giving the results for the 
(bulk) zero-temperature magnetization per site, magnetic susceptibility, density and compressibility: 



H 1 \ / ln(2) /8(/2 + 2a)\ 

m,uik = (l + ^y 3^(3) ) + 



1 / L\^{2) 8(/2 + 2a) 



27r2 V 21n(F) / \^ vr y 3C(3) ) 

21n(2) H"^ 1 \ 1 

+ ^ 4^ I " ln(//)j 7C(3)(/i + 2a) 

= ^TV 6C(3)i/2 + 2aj+-" • ^'-''^ 




The physical implications of ( 4.13| ) have been discussed extensively in the literature |43, 44 1. We 



merely note the divergences in the susceptibilities as we approach half-filling. 
4.1 Ground State Properties of the Impurity 

The impurity contribution to the ground state energy is of the same form as the surface energy of 



an open t-J chain with boundary chemical potential studied in [46| (more precisely the impurity 
contribution is of the same form as the boundary chemical potential dependent part of the surface 
energy). The parameter a plays a role similar to the boundary chemical potential in the open chain. 
The leading order impurity contributions for small bulk magnetic fields close to half-filling can thus 
be calculated in the same way as in |46| with the result 

+ ifa»i? 
' 9 \ 4.14 
1-||)+... ifa<B ^ ' 



Together with ( [4.q ) this gives the leading impurity contribution to the ground state energy. Note 
that by differentiating the impurity contribution to the ground state energy w.r.t. a it is possible to 
evaluate the expectation values of various operators at the impurity. However due to the complicated 
structure of these operators it is difficult to extract useful information from the expectation values 
and we therefore do not present these calculations here. 

For densities below and slightly above ric analytical results for particle number and magnetization 
can be obtained as well, whereas for generic values of a the integral equations can be solved only nu- 
merically. Below we first present the aforementioned analytical results and then turn to the numerical 
solution of the integral equations. 



4.1.1 Large a close to half-filling 

By taking the appropriate derivatives of the impurity contribution to the ground state energy we can 
evaluate the impurity contribution to the magnetization, particle number and susceptibilities. On 
physical grounds it is reasonable to assume that the impurity contributions to magnetization and 
particle number are concentrated in the vicinity of the impurity (see also Appendix B). We find 



M„ = _ I-—— + 



XH,a 




vrSay 3C(3) w->a ln(i7) V6C(3)(/Z + 2a) 
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djjL 



4 /2(/x + 2a) 
toV 3C(3) 



vra ^6C(3)[/i + 2aJ 



+ . . 

+ ... 



(4.15) 



By inspection of ( 4.15| ) we find that close to half- filling the impurity is on average almost doubly 
occupied and thus only weakly magnetized. The susceptibilities exhibit the same types of singularities 
as the bulk. 

In comparison to the bulk the following ratios are found: 



^ _A_ 2{j2 + 2a) XH,a 
rribulk TraV 3C(3) ' XH,bulk 



4 

7ra 



2{fjL+2a) 
3C(3) 



if -ff^ < A 



2 

aln(2) 



i-filling 



(4.16) 



The analytic (4.15) and numerical (Fig. ^ results show that the impurity is doubly occupied in 
the limit q — > oo. Thus the effective ground state hamiltonian in this case is a t-J model of L sites 
with twisted boundary conditions (see ( |2.12|) ) as represented by Fig. ^. This is consistent with the 
twisted boundary conditions occuring in ( p. 3D in this limit. 



4.1.2 Small q close to half-filling 

For small values of a such that a <^ B we find the following results for the impurity contributions to 
magnetization, particle number and susceptibilities 



XH,a 



H 
2^ 

1 

2^ 



1 

1 

a 



1 



2\n{H) 
1 

2\n{H) 



a 


)(■• 
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^ 4^4 





21n(2) 2{fi + 2a) 



3C(3) 
a 



+ ... 
1 



ln{H) J 2-kB {Jx + 2a) 



+ 



^^_^,21n(2) 



a ^21n(2) 2{jl + 2a) 
~V 3C(3) 



+ 



a 



ttB' tt v/6C(3)[A + 2a] 2TiB{ji + 2a) 



+ ... 



(4.17) 



Comparing this to (4.13) we see that the impurity behaves almost like an "ordinary" site of the lattice 
and the model reduces to a t- J chain on L -|- 1 sites in the limit a ^ (see ( 2.12| )) as represented by 
Fig. |. 



4.1.3 Densities below the critical electron density Uc 

Below the critical electron-density (|4.7] ) the particle number at the impurity-site is given by (see 
Fig. I) 



, , 2 /cot(^ 

Naine) = 1 arctan — 

IT \ 1 -|- a 



(4.18) 



In this regime only the states |0) and | |) are realized in the ground state, hence the hamiltonian 
7^l^2,3 simplifies to 



n 



1,2,3 



1 + 



-X. 



TT 



1 + - 



a 



^ [xl'x',^ + h.c.) + + +X]°X3°^ + h.c. 



(4.19) 
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4.1.4 Densities slightly above the critical electron density Uc 



Slightly above ric the impurity contributions to magnetization and particle number are found to be 

= 1 - - arctan — ^ + 

TT \l + aj \ ttBcJ ttBc 

11 f 2Ac \ f I \ a 

Ma = arctan IH h • • • (4.20) 

2 TT Vl + a/V T^BcJ 27rBc ^ ^ 



Combining ( |420| ) and (plO| ) we obtain the following result 

dNa. _ a + 1 - Na{nc) 

dug '^<=="c 1 — nc 

4.1.5 Numerical Results for general band filling and magnetic field 



(4.21) 



For general magnetic fields and band fillings the impurity magnetization and particle number can be 
determined by numerically solving the relevant integral equations. 

Some results for the impurity magnetization as a function of the magnetic field for different values 
of a and band fillings are shown in Fig. |5| and Fig. |6[ 

Only magnetic fields below the critical values are considered. For comparison the analytical results 
obtained above are included in the figures. The shape of the magnetization curves is entirely different 
from the ones obtained in the Kondo model where (as a function of H) there is a crossover from a 
linear behaviour for small fields to a constant value at large fields. 

From ( [4.18| ) one finds that the impurity magnetization at the critical magnetic field is an 
increasing function of a. Numerical results show that this behaviour persists for H < He, small a 
and not too large densities (see Fig. |5| and Fig. 0) . On the other hand the impurity is almost doubly 
occupied and hence only weakly magnetized for large a and magnetic fields sufficiently smaller than 
He- This leads to a decrease of the magnetization with a, which for small densities and small magnetic 
field fields is sh own i n Fig . ^ and Fig. 0. Near half-filling the magnetization is a decreasing function 
of a (see ( [4.15| ) and ( |4.17D ) for all magnetic fields not too close to He (see Fig. ^ and Fig. |^. 

The impurity particle number as a function of band filling for various values of a is shown in 
Fig. ^. For small a Na at first increases linearly with band-filling and then curves upwards, reaching 
2 at half-filling. For very small a the deviation from the = ng-behaviour occurs very close to 
half-filling. For large values of a and densities above ric increases rapidly with band-filling, then 
evens out and eventually reaches Na = 2. The crossover between these two regimes occurs roughly 
for a ~ 1 for small magnetic field . The magnetic-field dependence of the impurity particle number is 
shown in Fig. ^, where Na is shown for the case a = 5 and several values of the magnetic field. We 
see that apart from the shift in critical density ric the behaviour of Na is qualitatively unchanged as 
H is increased. The derivative of Na at ric is given by (4.21). 



4.2 Fermi velocities 

As we will now show the impurity contribution to the susceptibilities are related to a modification of 
the Fermi velocities by the impurity. The bulk Fermi velocities are defined by 



e'c{B) ^. e'M) 



" 27r pc{B) ' ' 27Tps{A) ' 
where the root densities Ps(A) and Pc(A) are solutions of the integral equations 

rA rB 

PsiX) = ai{X) - dp a2{X- p) ps{p) + I dp ai{\ - p) pc{p) 



(4.22) 



rA 

Pc{X) = / dpai{X- p) ps{p) . (4.23) 



A 
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The Fermi velocities can be calculated for small magnetic field near half filling (using the same tech- 
niques as above) leading to the results 



3C(3) 



B + 



V. = TT [1 



21n(2)S 



1 



+ 



21n(2)" ' " V vr (2\n{H/Ho)Y 

where Hq = ^/Stt^Jc. Near the critical electron density Uc we obtain 

_ 2 4Ac _ 8Ac 

~ ~ (1 + 4^2) arctan(2ylc)-Bc ^ "■' ~ 1 + 4^2 + - 

The respective Fermi velocities of the entire system (bulk+impurity) Vc and Vs are given by 

4(^/3) _ 



27r ^pp(Ap) + y^^'\Ap) 



P = S,C 



where Ag = A, Ac = B and satisfy the following integral equations 



A pB 

dp a2{\ - /x) p^Piii) + / dp ai(A - p) p^P{im) 

-A J-B 



p'^J'Hf^) = aa{n)+ I dp' ai{p - p) p^P {p) 



For later convenience we define the ratios 



f/3 



In the bulk the Fermi velocities are related to the susceptibilities via |44| 



1 



Xc,bulk 



vrn^ 
1 



XH,bulk . 

47r 



(Zee ~ '^Zsc)'^ ^ {Zcs — 2Zss)2 



Here Za/3 are the elements of the so-called dressed charge matrix (see e.g. 



Zee Zge 
Zes Zgg 



where ^^/^(A) are solutions of the integral equations 

















V 









j'^j^dp ai{\- p) - J^^ dpa2{\- p) 



iee{B) U{B) 

ies{A) iss{A) 



La dp ai{^- p) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



/ iceip) ise{p) \ 
V ics{p) Css{p) J 

Replacing the velocities in ( [4.29| ) by Ve and Vg we see that the impurity contribution to the magnetic 
susceptibility is given by 



XH,c 



1 

4^ 



{Zee ~ '2Zge)'^fe iZeg — 2Zgg)'^ f. 



+ 



(4.32) 



10 



There are two impurity contributions to the compressibihty: one is due to the change of the electron 
density rig, the other to the modification of the Fermi velocities 



Xc,a 



2Na 

DbulkL 



Xc,bulk ~l~ 



1 



imiL 



+ 



(4.33) 



The second part can be identified with -J—Q^ 



Ln'l d/j. 

Let us consider the two cases for which we presented analytical results above in more detail: 

• A ^ 1, B <^ 1 corresponding to a small bulk magnetic field and an almost half-filled band. The 
leading contributions to the dressed charge matrix are given by 



2 ln(2) 

Zr.r = l + —B 



1 - 



1 



41n(i7) - 21n(87r3/e; 



1 + 



1 



V2 V 4=HH/Ho] 



(4.34) 



where A,B,a are given by (4.12). With aid of the Wiener-Hopf technique we can obtain the 
ratios fc and fs as well. For a ^ B they are given by 



fc 



2 2a 

+ 



ln(2)a W(2)a 



+ •••, fs 



AB 



+ ... . 



(4.35) 



Inserting these results in ([4.32| ) and ( 4.33| ) we reproduce the results ( 4.15| ) (to ^ii^jj^))- For 
a <C -B we find 



/c = l + 



a 



21n(2)B2 y in(2) 



1 



+ 



a 



ttB 



(4.36) 



which leads to the same resluts as (|4.17| ). 

< ne — ric <C 1 corresponding to densities slightly above the critical density ric- We find 



Zr.. 



1+1 



1 



ttBcJ ttB^ 



+ 



Ar_ 



irBr 



Zss = l + ... 



+ .... 



(4.37) 



The ratios //? are given by 



2 arctan (|^) + vra ai+„(^e) 



2 arctan(2^c 



ai{Ac 



(4.38) 



5 The Impurity at Finite Temperatures 

The thermodynamics of the supersymmetric t-J model was studied by Schlottmann in The TBA 
equations for the dressed energies are the same in the presence of the impurity in complete analogy 
with e.g. so that we can simply quote the result from |3C] 

TT OO 
11 ^ X — ^ 

e = 27rai - // - — + rai *ln(l + e"T) -T2^ On * ln(l + e"~) 



n=l 



n 



2^1 + 27ra2 + Tai * ln(l + e"^) + Ta2 * ln(l + e"r) 

OO 

H - rin(l + e"^) + Tan * ln(l + e'^) + T J2 ^nm * ln(l + e 



(5.1) 



m=l 
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where 



(kve 



-iuiX 



coth(|||)(e 



-If (m~n)| 



f (m+n) 



The bulk free energy is given by |30| 

i^buik = -^'(0) - 2n, 
whereas the impurity contribution to the free energy can be cast in the form 

2a 



imp 



-2/i 



a + 2 



dX aa(A) 



^'(A) + rin(l + e" 



*{A) 
T 



(5.2) 



(5.3) 



(5.4) 



We note that in the zero temperature hmit this reproduces correctly the ground state energy ( [4 .2] ) . In 
the high-temperature limit the TBA equations (5.1) turn into algebraic equations that can be solved 
by Takahashi's method |48]. The leading terms of the high-temperature expansion are given by 



-^bulk 

L 



-Tin 
2a 



M H 

1 + eT2cosh — 



imp 



a + 2 



Tin 



1 + eT -\- eT 2, cosh — 
2T 



(5.5) 



We see that these yield the correct values of the entropy of a system of L sites with 3 degrees of freedom 
and one site with four degrees of freedom in the limit T ^ oo. We also note that the parameter a 
enters only in a trivial way into the leading term of the high-temperature expansion. By taking the 
appropriate derivatives we can compute the mean values of particle number and magnetization 



{D 



bulk/ 



2 cosh ^ 



imp/ 



exp(— ^) + 2 cosh ^ 
exp(^) + cosh ^ 
cosh(^) + cosh ^ ' 



("Ibulk) 



sinh ^ 



exp(— ^) + 2 cosh ^ 



(Mi 



sinh ^ 



imp/ 



cosh(^) + cosh ^ 



(5.6) 



Half-filling corresponds to the limit // — > oo, in which the impurity is on average doubly occupied and 
unmagnetized whereas the bulk exhibits a magnetization per site of |tanh^. 



6 Low Temperature Specific Heat 

In order to calculate the contributions of the impurity to the low-temperature specific heat we need to 
consider different characteristic regions (see Fig. 10) of the t-J model as it was done for the Hubbard 



model by Takahashi As in we assume H ^ T so we can neglect the effects of the string 

solutions. (An alternative approach which overcomes this restriction has recently been used in |45] to 
compute bulk thermodynamic properties in the t-J model). 

Region A: The region is characterized by fx < — -y- For T = the electron density is zero as 
A = B = oo. The low temperature free energy per site is given by 

AFbuik _ Fbuik - i^buik(T = 0) _ r ,,,,,,, ^ 



L L 



V2^ 2^ 
The impurity contribution to the low temperature free energy is 



T dAa2(A)ln(l + e— )-r / dA ai(A) ln(l + ) 

J —oo J —oo 

2T f°° I , 2m T 1 , ,^ 

/ dX — ln(l + er e ) / dX —r ln(l + e t e a2t ) 

vr Jo A^ vr Jo X^ 

1 ^•^/O 2m 1 „q/o M+g/2 

T^/^e- —T^'^e^r^ (6.1) 



dX a,(A) ln(l + e— ) « --^T^/^eT . (6.2) 
-oo 2v ivr 
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Region B: For T = this region, characterized by —H/2 < fi < H/2, is the ferromagnetic phase 
with electron density varrying as < ng < 1. The right boundary hne is defined by i? = oo the left 
one by ^ = -B = oo. The low temperature free energies are given by 



bulk 



L 



ttT' 1 
6 



6 V, 



(6.3) 



Region C: In this region, characterized hy < A,B < oo, the electron density varries between 
< rie < 1 at zero temperature. The right boundary line defined hy B = can be calculated with 



aid of ( [4.12| ) for small magnetic field and with an iterative solution of 
temperature free energies are given by 



bulk 



6 



1 1 



6 



fs ^ fc 

Vs Vc. 



for ~ 4. The low 



(6.4) 



Region D: This region is characterized by ^ > H/2 > 2. For T = we obtain the ferromagnetic 
half-filled band. The low temperature free energy is given by 



AR 



bulk 



L 



T dX ai(A)ln(l + e" 



4T r°° , 4-H. -16a2 . 

/ d\ ln(l + e T e T 

vr Jo 



The impurity contribution to the low temperature free energy is 



AF„ 



r / d\ aQ(A)ln(l + 



HI2-U. 



--l=r3/2e^ . (6.5) 



(6.6) 



Region E: For T = this region is half-filled and non ferromagnetic. The low temperature free 
energy of the bulk is given by 

The impurity contribution in this region can not be calculated in closed form. In the two limiting 
cases A <^ 1 and ^ ^ 1 we obtain the following results 



2 ln(2)-)j+2a 

=e T 



AF„ 



'_rp3/2 



a 



A/37rC(3)/2-7r3a 

H/2-M+^(4-g)=^/^ 



for ^ > 1 
for ^ < 1 



(6.8) 



Wilson ratio in Region C: The specific heat of the bulk and the impurity in region C are given by 



vrT / 1 



CvJyulk — 'TT ' 



1 



^v.a — TT^ 1 



3L \v 



(6.9) 



In deriving these results we assumed that H ^ T. Defining 

XH,a ,bulk 



R 



-2Zs. 
-2Zs 



( Zee ^Zsc \ 
Zcs 2Zss J 



fs Vs 




1+ ^ 


_|_ Vc_ 

Vs 




Ic , Vc 
Js Vs_ 



(6.10) 



and if we assume that the limits T — > and H ^ commute Q we can calculate a "Wilson ratio" 

1 + ^ 



Rw 



lim R 



fc _l_ Vc 
fs Vs 



(6.11) 



'*This holds in all known cases for the specific heat in integrable spin chains where the ground state contains only real 
roots of the Bethe equations. The assumption is also supported by the findings of where it was shown to be true for 
the bulk specific heat. 
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Unlike for the case of the Kondo model (where spin and charge degrees of freedom decouple and the 
impurity couples only to the spin) there is no reason to expect Rw to be universal for the present 
model. 

From our analytical calculations we find for B <^ 1, B <^ a that Rw = ^ — <^ 1. Near the 
empty band the ratio Rw tends to one because Vc approaches zero. Numerical calculations show that 
the limiting value one is most rapidly approached for small values of a (see Fig. O). For comparison 



we quote the result for a Kondo impurity in a Luttinger liquid p|, for which Rw = ^{1 + 



7 Transport Equations 

Following Shastry and Sutherland [^] (see also [^]) we will now calculate spin- and charge stiffnesses 
from the finite-size corrections to the ground state energy of the model with twisted boundary con- 
ditions. For the t-J model the bulk stiffnesses were determined in |4^. Our analysis follows closely 
the discussion given in |44|. Our goal is to evaluate the ground state energy as a function of the twist 
angles (p-^ and Imposing twisted boundary conditions the BAE ( |3.3| ) are modified in the following 
way 

11 , (1) , i ' -J 




'y=~\_ /\'y 

11 (1) i 11 (1) (1) ' 7 ^'•••'^n • K'-^) 

j = l Af — Aj — 2 (3=1 A-f — — ^ 

For technical reasons it is convenient to use a different representation of the BAE first introduced 
by Sutherland for the t-J model without impurity |29|. This can be obtained by a particle-hole 
transformation in the space of rapidities which is done in Appendix A. The final BAE are given by |^ 

L 




A/'i)+Wi-l A/(i) , ,(1) , 



X-j — Xk — i A — A^^-' 

k=l ■' " 7 = 1 '^J ^7 



M'i)+Wi-1 (1) , (1 



n a)';'? - ^"^'^T?T)^ 7 = l,...,M«=i + 2-iV,-^,, (7.2) 

3 = 1 ^1 + 2 -^7 * 2 



A«'' = AKi^-i , i = l...n, (7.3) 



where the twist-angles are given by (ps = (t>'[ — 4>i and (pc = (pi |44]. 

The equations ( [7. 21 ) can be simplified by making use of the 'string-hypothesis'|^, which states that 

for L ^ oo all solutions are composed of real A^-^^^'s whereas the A's are distributed in the complex 
plane according to the description 

'n + l 

where (3 = 1... Mn labels different 'strings' of length n and A^ is real. The imaginary parts of the A's 
can now be eliminated from ( [7.2p via ( [7.31 ). Taking the logarithm of the resulting equations (for M„ 
strings (|7j) of length n and M^) A^^^'s (note that nM^ = L -M - iV^)) we arrive at 

T^." = «(|)4E«-(^5-A")+ji:'<'(^^i^)4,/3=i...M,. 

(m7) ')'=1 

^J, = yE^K^^-A^) + >(4^) + T'7 = l...M«, (7.4) 



Considering the half-filled case A^^ = we see that the impurity leaves the BAE (7.2) unchanged, which is not 
immediately obvious from the other set of BAE (7.1). 

®Note that we do not consider string solutions to the BAE in order to determine the stiffnesses and the results are 
thus independent of the precise form of the strings. 
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where and J-y are integer or half-odd integer numbers, 9{x) = 2arctan(2x) and 

en,U^) = (1 - + 2 o{- — ^— ) + . . . + 2 e{ / ) + e{^—) . (7.5) 

|n — m| |n — ml+z n + m — 2 n + m 

For vanishing twist angles the ranges of the "quantum numbers" and are given by 

oo 

< -(L + M« + A/„-2 ^ min{m,n}M„- 1) , 

m=l 

1 oo 

< -iY^Mn-l). (7.6) 

^ n=l 

Ground state and excitations can now be constructed by specifying sets of integer (half-odd inte- 
ger) numbers IJ^ and and turning equations ( |7.4D into sets of coupled integral equations in the 
thermodynamic limit. The antiferromagnetic ground state for zero twist angles is obtained by filling 
consecutive quantum numbers ijj and symmetrically around zero, which corresponds to filling two 
Fermi seas of spin and charge degrees of freedom respectively. The effect of an infinitesimally small 
flux is to shift the distribution of roots (i.e. the rapidities in the Fermi seas) by a constant amount. 
This shift leads to a twist-angle dependent contribution to the ground state energy. The ground state 
(in the presence of flux) is described in terms of the root densities /5o(A), which are solutions of the 
integral equations 

r'At rA+ 
PsW = ai{X) - I dfi a2{X - PsifJ-) + I d^x ai{\ - ^j) p^^) + 0{L ) 
JAj Ja- 

Pcip) = J^! dp' a,{p - p') psip') + ^ + 0{L-^). (7.7) 

Here are the Fermi points for the finite system in the presence of the flux. We denote the Fermi 
points for the infinite system without flux by ±A^_o (note that the distribution of roots is symmetric 
around zero in that case). For further convenience we define the quantities 

Dp = -o /--/ dp pfsip) . (7.8) 



2Ja 

In order to evaluate the stiffnesses we need to consider infinitesimal flux only, which yields a correction 
of order to the ground state energy. Following through the standard steps |5^ and taking into 
account the -^-corrections from the infinitesimal flux we obtain^] 

Ei<Ps, ct>c, a) = Le(A^,o) + /(A^o) - ^{vs + Vc) + I-kLD^T^VT^D + o(-^), (7.9) 

where Z is the dressed charge matrix ( [4.30 ), D = {Dc,Ds)~^ and V = diag(uc,fs). Here e(A^o) and 
/(A^ o) S'l's the ground state energy per site and impurity energy in the infinite system without flux. 
The quantity -D/3 is given by 

= 1^1 ■ 

In order to study the charge and spin conductivities and the respective currents we consider the 
energy-difference AE = E{(ps, 4>c, ct) — E{0, 0, a). It can be cast in the form 

AE = Y^oiE>ai34>i3 + subleading terms. (7.11) 



'^Similar expressions are obtained for the corrections to excited state energies. 
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According to [50| the charge (spin) stiffness D^''^ {D^'^^) is defined as 



^^'^ = = ^As = 0, a) , I)M = ^^AE{cP, = <t>As = -20, a). (7.12) 



Using the expressions (|7l0| ) in (|7j) and ( fml) we find that the stiff'nesses are not modified by the 
impurity to leading order in . It is clear from the above analysis that there are corrections due to the 
impurity in the subleading terms. The precise form of these expressions is not of particular interest 
from a physical point of view and as the extension of the above finite-size analysis to the subleading 
orders is rather difficult we refrain from determining them. The important point is that despite the 
presence of the impurity the stiffnesses are still finite and the dc-conductivity is thus infinite. We 
believe that this fact is due to the integrability and the related absence of backscattering. 

This means that the integrable impurity considered here is of a completely different nature than 
the "weak link"-type potential impurity considered in [|l|, ||, |3|, as the electron-electron interactions 
are repulsive in the supersymmetric t-J model a weak link drives the system to a strong coupling fixed 
point characterized by the vanishing of the conductivity. 

7.1 Stiffnesses for finite density of impurities 

Let us now turn to transport properties for the system with a finite density of impurities. The necessary 
steps are the same as above, the main difference being the change of integral equations describing the 
ground state from (|7.7|) to 



r-A+ M + 

Ps{X) = (1 - ni)ai(A) - / dn a2{X - PsifJ-) + 1 d^i ai[\ - fi) pcip) , 



a: 



Pcilj) = niaa{p)+ I dn' aiiiJ. - fi') ps{p') , (7.13) 
Ja, 

where is the concentration of impurities. As the integral equations for the dressed energies remain 
unchanged, the dressed charge is not modified and A{p,H) and B(p,H) are not changed. However, 
the presence of a finite density of impurities leads to changes in the electron density, which is now 
given by ne = (1 — nj)-D;,u;fc -|- UiNa and the Fermi velocities, which are found to be of the form 

1 + - 1) 



Here Vj^ are the velocities of the normal t-J-chain and the fp are defined in (|]2|). From ( pl|) the 
following form of the stiffnesses is easily deduced 

D^f^ = ^ {vcZl + VsZi) , = ^ {vc{Z,, - 2Zscf + UZcs - 2Zssf) . (7.15) 

Using the results of the above sections we can evaluate these expressions analytically for small magnetic 
field close to "maximal filling" (t-J sites singly occupied, impurities doubly occupied) and near the 
critical electron density ric- This is done in the following two subsections. Finally we present numerical 
results for the general case. 

7.1.1 Stiffnesses for small magnetic field near maximal filling 

Close to maximal filling and for a ^ B the leading term of the charge stiffness is found to be 

~ ^2 _ 3C(3) B 

^ ^ -h \^i^2)^ + ln^(2)a '-J 
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Combining (|7.16| ) with the result for the electron density we obtain the following limiting value for 
the slope of the charge stiffness as a function of density 



drip 



|ne=l+ni 



3C(3) 



81n2(2) 



1 + n, 



Vln(2)a 



_l_ 2a 



In the strong-coupling limit this turns into 



lim 



dD^p"! 



+00 Qup 



\ne=l+ni 



3C(3) 



1 dD 



ip) 
t.j 



;ln2(2)fl-n,V 



|ne=l 



The leading term of the spin stiffness for B -^1 and q ^ i? is given by 

1 V 1 



.2 Vs 



2-K 

In the limit a — > cx) we find 



vr V 4ln{H/Ho)J l + ni(/,-l 



dup 



|ne=l+ni = vr 1 



1 dD 



tj 



2ln{H/Ho)J {l-my (1 



n,; 



n,=l 



(7.17) 



(7.18) 



(7.19) 



(7.20) 



7.1.2 Stiffnesses slightly above Uc 

As pointed out above it is possible to derive analytic expressions for the stiffnesses for densities slightly 
above ric- However the resulting expressions are found to be rather complicated so that we refrain 
from listing them here. 

The derivative of the spin-stiffness with respect to the electron density at rip = is always 
positive, taking its maximum at a = and its minimum in the limit a ^ 00. 

The derivative of the charge-stiffness with respect to the electron density at rip = changes sign 
as a function of a (see Fig. |l^) . 



7.1.3 Numerical Results 

The results for the charge stiffness in systems with impurity densities rii = 0.2 and two different values 
of the bulk field H are depicted in Fig. 13 and Fig. The charge stiffness for = 0.5 is shown 
in Fig. 15. For comparison we plot the result for the charge-stiffness D^j for the t-J model without 



impurities as calculated in |44]. We note that the maximal allowed band- filling is larger than one as the 
impurity sites can be doubly occupied. We see that for small band-fillings above the critical density 
the charge-stiffness is reduced as compared to the pure t-J case. For larger band fillings the stiffness 
is found to increase in the presence of impurities. This is easily understood: due to the constraint of 
single occupancy the stiffness vanishes as we approach half-filling in the t-J chain. The impurity sites 
can be doubly occupied, which gives the electrons "space to move" and leads to an increase in the 
stiffness. For large fillings the stiffness increases with increasing a because (as can be deduced from 
the a ^ 00 limit) the impurity sites become (on average) closer and closer to being doubly occupied, 
which again makes it easier for the electrons to move along the t-J sites. Last but not least let us 
discuss the limiting case a = 00 at impurity density rii, i.e. there are Lui impurity sites and L{1 — rii) 
t-J sites. For electron densities He smaller than rii the (spin-up) electrons (in the ground state) occupy 
only impurity sites which do not interact with the t-J sites. Thus the stiffness is identically zero. For 
electron densities < ng < + (1 — ni)nc the saturated ferromagnetic ground state on the t-J sites 
is formed whereas all impurity sites are singly occupied. The stiffness is completely determined by 
the t-J sites. For densities in the interval nj + (1 — nj)nc < rig < 2nj + (1 — ni)nc the impurity sites 
become doubly occupied and the stiffness does not change. For rig > 2nj + (1 — ni)nc all impurity 
sites are doubly occupied and the stiffness follows (up to a rescaling by the t-J curve above the 

critical density ric (see ( [7.20] ) and ( |7.1^ 



The spin-stiffness for impurity density 0.2 is shown in Fig. 16. We see that the stiffness is decreased 
at low fillings (this decrease is more pronounced for larger values of a) and approaches the "pure" t-J 
value for large fillings. The behaviour in the limit a — > 00 is the same as for the charge stiffness. 
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8 Impurity Phase-Shifts 



In this section we evaluate the phase shifts acquired by the elementary excitations, holons and spinous, 
when scattering off the impurity. The results give a good measure of the effects of the impurity on 
excited states. In particular we can infer from the phase-shifts how the impurity couples to spin and 
charge degrees of freedom. We start by briefly reviewing some important facts about the low-lying 



excitations in the t-J model [31]. The elementary excitations are collective modes of spin or charge 
degrees of freedom. The spin excitations are called spinous and carry spin it^ and zero electric charge. 
They are very similar to the spin-waves in the Heisenberg XXX chain. The charge excitations are 
called holons and antiholons, carry zero spin and charge ^e. Thus holons correspond to physical holes. 
At half-filling only holons can be excited as the charge Fermi sea is completely empty. The excitation 



energies are given by ec,s defined in (4.3). The respective physical momenta are given in terms of the 



solutions of the following set of coupled integral equations 

PsW = -^(A) - / du a2{\ - v)^s{v) + I dv ai{\ - u) , 
J-A J-B 

rA 

p,(A) = / dz. ai(A-i^) p,(i/). (8.1) 

J-A 

The mometum of e.g. a holon-antiholon excitation is given by Pec = Pc{A.P) - Pc(A'') where A^ and 
A'^ are the spectral parameters of the holon and antiholon respectively. We thus would define the 
physical holon momentum as pc(A^) = pc(AP) — pc(i?). At half-filling the spinon (ps) and holon (pc) 
momenta are given by 

Ps(A) = 2 arctan(exp(7rA)) — TT , 

PciX) = ^ + zlnf^^^^^i±|lV (8.2) 

I'-y J 2 ^ I r(i±^)r(i-f)y ^ ^ 



The scattering matrix has been determined by means of Korepin's method [39, p3] in [31|. At 
half-filling the spinon-spinon S-matrix S{X), the spinon-holon (sc) and holon-holon (cc) scattering 
phases are given by 

.(A) ^ ,n^-^)^il^%)f A 



r(i^)r(i-f) Va 

l + ie^^ r(i+iA) r(l - — ) 

exp(i(/p,,(A)) = , eMi^ccW) = r(i-}A)r(i I \) ^^-^^ 

where / and P are the 4x4 identity and permutation matrices respectively. Below half-filling the 
S-matrices are given in terms of the solution of integral equations. 

The impurity phase-shifts can be computed by the standard method of Korepin [^], Andrei et. al. 



55 [. In the most general case of a bulk magnetic field and arbitrary filling factor the phase-shifts 
can be expressed only in terms of the solution of a set of coupled integral equations, the analysis of 
which is rather difficult. We therefore constrain ourselves to the case of a microscopic number of holes 
in the half-filled ground state in the absence of a bulk magnetic field. 

The basic ingredient for computing impurity phase-shifts is the quantization condition for fac- 
torized scattering of two particles with rapidities Ai^2 on a ring of length N (including the impurity 
site) 

exp(iAfp(Ai))5(Ai - \2)e'''^^^^^ = 1 , (8.4) 

where p{\) is the expression for the physical momentum of the corresponding (infinite) periodic system, 
iS'(A) are the bulk scattering matrices for scattering of particles 1 and 2, and is the phase-shift 

acquired by particle 1 when scattering off the impurity. We note that the condition ( ^.4D incorpo- 
rates the fact that there is no backscattering at the impurity. For the present model the absence of 
backscattering follows from the conservation laws for the rapidities: although momentum is not a good 
quantum number for the ring with impurity, excited states can still be characterized by the rapidity 
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variables (see below for an example). We would expect that if the impurity contained a backscattering 
term mixing of states with different rapidities would occur. This is not the case in the present model 
which indicates the absence of backscatterin g. W e note that the absence of backscattering on the 
level of the "bare" Bethe Ansatz equations ( |3.3|) (which describe the scattering of excitations over 
the empty ground state) is not sufficient to deduce the absence of backscattering over the antiferro- 
magnetic ground state because the impurity gets dressed by the holons and spinous in the ground 
state Fermi seas, and the two-particle scattering processes between holons and spinous do contain 
backscattering terms. Therefore the treatment of |^6[ does not apply in the present case. 

In what follows we will extract the holon and spinon impurity phase-shifts from the spinon-holon 
scattering state, for which the condition (8.4) turns into scalar equations for the scattering phases, 
which after taking the logarithm read 



Np,{X'') + ipsciX'' - A*') +M>^'") = mod 2tt 
Npc{AP) + ipsci^P - \^) + Vc(A^) = mod 27r 



.5) 



Here A'' and A^ are the rapidities of the spinon and holon respectively. Comparing these conditions 
with certain quantities ("counting functions") that can be calculated from the Bethe Ansatz solution 
one can then read off the boundary phase-shifts ips,c- 

Let us start by constructing the half-filled antiferromagnetic ground state for even length L of the 
host chain, where we furthermore assume that is even as well. The ground state is obtained by 



choosing Mi 



in ( [7.4| ) and filling the half-odd integers l\ symmetrically around zero. 



In the thermodynamic limit this corresponds to filling a Fermi sea of rapidities \\ between — cxd and 



oo, where the root density Ps(Al) = lim^,. 



(L+l)(Ai^^-Ai) 



is given in terms of the integral equation 



p,(A) = ai(A) 1 



1 



L + 1 



du a2(A - v) ps 



16) 



The spinon-holon scattering state characterized by choosing Mi 



1 in the Bethe 



equations (7.4). There are ■§ + 1 vacancies for the integers and thus one hole in the Fermi sea of 



A^. We denote the rapidity corresponding to this hole by A^. The rapidity corresponding to the holon 
is denoted by A^. The Bethe equations read 



27r 

L + l 



J 



L + 



A/3 — A/3/ 



1 



/3'=1 

1 „,AP. 



L + 



/3=1 



+ 



L + 1 



,A^-A^ 



L + l 



e{KP - A^), 



) + e{\p - An 



(8.7) 



where J is a half-odd integer number. In the limit L ^ oo the distribution of roots A/3 is described 
by a single integral equation for the density of roots Ps(A), which is of Wiener-Hopf form but cannot 
be solved in a form sufficiently explicit for the purpose of determining the impurity phase-shifts. The 
main complication is that we need to take into account all contributions of order and thus must 
deal with the fact that the roots are distributed not between —00 and 00 but between two finite, 
-L-dependent values —A and A. It can however be checked numerically that making the assumption 
that the contributions due to the shift of integration boundaries will be of higher order in as far 
as the impurity phase-shifts are concerned (and thus taking A = 00) yields the correct result. The 
integral equation then can be solved by Fourier transform 

fTsiu;) = Goiu;) + {Gi(u;)e-^' - Go{u;)[l - 26^^^"} , (8.8) 

where ^(w) is the Fourier transform of Ps(w) and where Gx{uj) = '^2cosM^^f • 
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For the further analysis it is convenient to define counting functions Zs{X) and Zc{X) 



+1 



27r 



/3=1 



1 

+ 2^ 



0{ 



) + e{x-An 



1 ^+ 1 r A 

-y 0(A-A;3)-:^ ^(A-A'') + 0(- 
zvr ZTT a 



/3=1 



(8.9) 



Note that for any root e.g. Aq, of ( |8.7D the counting function takes the integer value Zs{Xa) = la 
by construction. In the thermodynamic limit times the derivative of Zs yields the distribution 
function of rapidities PsiX). Straightforward integration of the density ps{X) yields the following results 
for the counting functions in the thermodynamic limit evaluated at the rapidities of the spinon and 
holon respectively 



27rz,(A'^) 
-27rze(AP) 



(L + l)ps(A^) + ifsciX'' - AP) + 0,(A'^) = mod 2tt , 
{L + l)vc{Kn + Vsci^" - A^) + 0c(AP) = vr mod 2tt , 



i.W) 



where ps,c{ X) are the spinon/holon momenta ( ^.21) , ipsc{X) is the bulk phase-shift for spinon-holon 
scattering (|8.3|) , and 

<l)sW = -Ps{X) , M^) = -Pc{X)+7T-9{-) . (8.11) 



a 



Prom these equatio ns w e can now infer the boundary phase shifts by comparing them with the quan- 
tization condition ( ^.51) , which yields 



c 



r(i + f)r(^ 



2 ^ 



-c 



-■-14 2 ' ^4 ' 2 

1 a-2iA r( 



4 2 ^ V4 I 2 ^ 

1 _ iA> 

2 2^ 



r(i + 



a + 2ix r(i + f)r(i 



2 ' 



^.12) 



where C is an overall constant factor of unit modulus that cannot be determined within the Bethe 
Ansatz framework. Setting C = —i we find that 



-ips(A) 



P«V'c{A) 



a 



2iX 



a + 2iA 



-'tpc(A) 



^.13) 



where ps and pc are the spinon and holon momenta at half-filling respectively. This result is inter- 
preted in the following way: for the half-filled band doped with a finite number of holes the impurity 
site essentially decouples from the host chain in the sense that spinous and holons bypass it, which 
effectively shortens the lattice by one site (see Fig. ^! 

For the spinous this is the complete picture, whereas the holons still acquire a phase shift due to 
the fact that the impurity site is charged (recall that it is on average almost doubly occupied) and 
therefore interacts with the holons passing it by. The holon scattering phase has a pole at A = i^, 
which for —2 < a < —1 lies on the physical sheet and therefore corresponds to an impurity bound 
state. The restriction a < — 1 is imposed in order to have a hermitean hamiltonian The impurity 
therefore has the interesting property to lead to a holon bound state for sufficiently small negative a 
at half-filling. 



9 Conclusion 

In this paper we have studied the effects of an integrable impurity in a periodic t-J model. The 
impurity couples to both spin and charge degrees of freedom and the coupling strength a can be 
varied continuously without losing integrability. The two limiting cases a = and a = oo have been 
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shown to be described by effective (ground state) hamiltonians of a t-J model with one extra site, and 
a decoupled impurity in a t-J chain with twisted boundary conditions. 

At zero temperature we have calculated the impurity magnetization and particle number for arbi- 
trary band filling and bulk magnetic field. The impurity susceptibilities have been shown to exhibit the 
same types of singularities as the corresponding bulk susceptibilities. Similarly the low-temperature 
specific heat of impurity and bulk have the same temperature dependence. Transport properties have 
been determined through the calculation of spin and charge stiffnesses and finally the impurity phase 
shifts have been calculated for the half-filled band. 

The supersymmetric t-J model belongs to the class of Luttinger liquids with repulsive electron- 
electron interactions. The effects of potential impurities of the "weak- link" type were first investigated 
in ||l|, |2|. It was found that the system fiows to a strong-coupling fixed point characterized by the 
vanishing of the dc conductivity. The physics of the impurity studied here is quite different: the 
dc-conductivity is unchanged by the presence of a single impurity. 

As argued above the type of impurity considered here does not seem to contain backscattering 
terms on the level of the dressed excitations (holons and spinous). It would be interesting to verify this 
assertion by explicitly constructing the continuum limit of the model. However, due to the complicated 
structure of the impurity hamiltonian this is a difficult undertaking. The argument given above 
suggest that it is impossible to construct an impurity model containing backscattering off the impurity 
by means of the Quantum Inverse Scattering Method through the standard intertwining relation 
"i?TT = TTR": the rapidities of the elementary excitations will always be conserved quantities and 
are not affected by the scattering off the impurity. Clearly "generic" impurities ought to contain 
backscattering as only special potentials are reflectionless. From that point of view the integrable 
impurity considered in the present work is very special. The situation is similar to the (multichannel) 
Kondo model (viewed as a 1 — d system). One may speculate that like for the case of a Kondo- impurity 
in a Luttinger liquid a backscattering term will drive the system to a new fixed point |jl2| so that the 
integrable impurity would represent an unstable fixed point in the sense of the renormalisation group. 
This is known to be the case for the spin system of Andrei and Johannesson |57]. As we have seen the 
integrable impurity nevertheless leads to interesting physical consequences. 
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A Transformation of the BAE 



To show the equivalence of the two sets of BAE ( |7.1| ) and (^^) we use a method due to Woynarovich 
| p§| ] and Bares et aJ||44| . We express the second set of ( |7.2| ) as a polynomial of degree 



n 



{w — \j ){w — i—) — e 



M(i)+Ar|-1 
13=1 



(A.l) 



and identify the first M^^'^ roots of ( AJ) wi, . . . , Wj^,j(i) with Aj^^^ , • • • , ■ Using the residue theorem 
we obtain the following equality: 



A/ — A^"^"* -I- 1 
A/ - A^-^^ - § 



y — 



\ 27ri /c 



dz— In 



A; - 2 + 
\v-z-% 



2 J 



-ln(P(.)) 



(A.2) 
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where Cj is a small contour including A^^^^ . Deforming the contour and denoting the Ni other roots of 



( \AJ\ ) with w[,... 



w 



we arrive at the following equality 



+2 



"A/ 




i ' 
2 


+ - In 

i 








i 

2. 





(A.3) 



where the last term comes from integration aorund the branch cut extending from z„ = \i + i/2 to 



A^ — i/2. Using the form of P{w) and substituting ( A.3| ) into the first equation of ( |7.2| ) we obtain 



the first equation of (7.1) with the according twist-angle 

L 




A. 



.A,+ 



-i{4'c+(t>s) 



n 



A,- 



1 Aj 



(A.4) 



The second equation of ( |7.1| ) can be obtained by the same steps starting with the first equation of 
( [7.2[ ). The twist-angles are related by (j)s = (^^ — 4>i ^"^^ 0c = (^i 0] shown that 

the BAE (7.2) for the t-J model without impurity can be obtained by means of the QISM starting 
with a fermionic vacuum with all spins up. The corresponding vacuum state of the impurity model is 
given by a bosonic doubly occupied impurity site and all other sites occupied with spin up electrons. 
The algebraic Bethe-Ansatz starting from this vacuum can also be constructed. 



B The three site model 

The Bethe ansatz states do not form the complete set of eigenstates of the system but are the highest- 
weight states of the gl{2\\) superalgebra (taking the Lai solutions of the BAE). Complementing the 
Bethe ansatz states with those obtained by the action of the gl{2\\) lowering operators one obtains 
additional eigenstates. The completeness of this extended Bethe ansatz has been proven for some 
models as the spin-^ Heisenberg chain, the supersymmetric t-J model and the Hubbard model 



60 1 . In this appendix we present a completeness analysis for impurity system considered here on a 



chain with three sites. This nontrivial example shows that the picture of |33, 59, 60 1 seems to hold in 
the present model as well. A detailed analysis of the general case is outside the scope of the present 
paper. 

We need to consider the action (on states given by the Bethe Ansatz) of the spin lowering operator 
= Y^i=l and the supersymmetry operators Q\ {Qa in the Sutherland case), which are given 

by 

L+l L+1 

Qi= + ./^^xf + ^xl\ gT= E x^^ + ^/^xl'^ -^xi\ (B.i) 

They are seen to satisfy the commutation relations 

{Qt,Qi} = 0, Q^ = 0, [H,Q,]=0. (B.2) 



Mil 



The BA states obtained by the Lai solution starting with empty sites are characterized by Q(j\^ La 
0. The respective Sutherland solutions by Q]j\'^Suth.) = 0. Solving the BAE (7J) and ( |7.2|) with 



vanishing twist angles for the simplest case L = 2 (recall that L is the length of the host chain) 
and then constructing the corresponding gl(2\l) multiplet by acting with all possible combinations of 

raising generators we obtain the following complete set of eigenstates (A = ^1/ §^) 
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Energy 


u 

tl 


T3 A T tin 

dA Lai 


BA Sutherland 





4 


vacuum 
l*o> = |0) 


Ai = -A2 = A Ai^^ = -A2^^ = ^ ^ 

\^o) = Q]Q\\^o) 


6+2a 
a+2 


8 


Ai = A 9 = —\ yja + 1 + i\/a + 2> 

l^i) =5*1 Too)+5|ooT) + |OTO) 


Ai = A A^^^ = A^ 
|«>i) =Q{qI|^i) 


&+2a 
a+2 


8 


Ai = -A 

|^2) = 5l TOO) +5*100 T) + |o TO) 


Ai = -A aP = -A^ 

1^2) = q{qIi^2) 


12+Aa 
a+2 


12 


Ai = — A2 = A 

1^3) = 1 TT 0) - Va + 1| T T) + |o TT) 


vacuum 

I^3)=Q{qI|^3) = IT2T) 


4 
a+2 


4 


Ai = — A2 = § A]^ = 

|^4)=QTQil^4) 


Ai = 

1^4) = 1 T 2 i) - li 2 T) 



The Sutherland solutions for the states |<&o,i.2) and the Lai solution for [^'4) Qdo not exist in the 

case Q = 0. In this case the Qa^ commute with X^^ and the states decompose into the 27 states 
corresponding to the Lai-states without double occupancy and the 9 states of the Sutherland solution 
with doubly occupied imp urity site 0. This is in agreement with the form of the hamiltonian in the 
limit a ^ given in (|2.12| ). 

Let us conclude this appendix with some simple considerations concerning the question of whether 
the impurity contributions to particle number and magnetization are indeed concentrated at the 
impurity. The lowest energy state with one electron on a lattice of arbitrary length is given by (5t|0). 
Using this fact we are able to directly compute the electron densi ty at the impurity for this state and 
we find that (71-2) = 2,+i°q, ■ order to compare this with ( 4.18 ) we need to take into account that 
( [4.18| ) is obtained in the thermodynamic limit, i.e. we need to take L ^ a + 1. We then find that 



Na = ^-T^t which is in agreement with (n2). 



Similarly the lowest energy state above the critical density ric is given by Q^^\^ Ferr) ■ By the action 
of the Q-operator a spin-down electron is generated with probability (a + 1) x (1 — No) and a doubly 
occupied impurity-site with probability a x ■ Taking into account the normalization-factor i^^-^^ ^ 
we obtain the following result 



AiVg 



LAiV„ 



1 



(B.3) 



This coincides with ( [4.21| ). From these simple examples we deduce that the assumption that the 
impurity contributions to magnetization and particle number are located at the impurity is a very 
reasonable one. 



C gl{2\l) Invariance of the Model 

In this appendix we show by explicit computation that the model ( |2.8[) is 5(/(2|l)-invariant. We start 
by expanding i?-matrix and L-operators around infinite spectral parameters 

fl33(A) ~ no° + ^(/-nr , 

£^3(A) ~ /0" + ^(n-/)0" , 

Cm{\) ~ /°2 + l(£-(2 + a)/r , (C.l) 



The solution of the BAE exists but the norm of the corresponding state vanishes. 
^The muhiplet of dimension 12 sphts in one 7 and one 5 dimensional multiplet as Q|(a = 0)| lit) = 
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where we denoted the auxihary space by 0, n labels the quantum spaces over the t-J-like sites, and 
the impurity sits at site 2. This leads to the following expansion of the monodromy matrix 



r(A) ~/ + ^ 

A 



Y^iU - ir + iC-{2 + a)I) 



02 



Inserting and (|C.2D into the intertwining relation 

ii33(A - /i) (T(A) ® T{fi)) = (r(/i) T(A)) i?33(A - fi) 
we obtain the following equations 



(C.2) 



(C.3) 



(C.4) 



Setting a = b', multiplying (|C.^ ) by (— and then summing over a we arrive at 

[y'"^',r(/x)]=0, (C.5) 
where r(/i) is the transfer matrix of the system. Dropping some constants we therefore find that 

[Qab,r{fi)]=0 , Qab = Y.^n a, 6 = 1 . . . 3, (C.6) 

n 

where we use the correspondences 1 <->t, 2 <->J,, 3 <-> for X""^ . The operators Qab form a complete 
set of generators for gl{2\l), which establishes the invariance of our model. 
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Figure Captions 
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Figure 1: The circles denote tlie normal t-J-sites. The square denotes the impurity placed at the 
second site. 
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Figure 2: Hamiltonian in the limit a oo. As in the half filled case the double occupied impurity 
site decouples of the bulk-chain which has to be solved with twisted boundary conditions. 
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Figure 3: Effective ground state hamiltonian in the limit a — 0. The impurity site is just an extra 
site of the t-J model. (This is correct for all particle densities below half filling where the impurity 
site becomes doubly occupied.) 
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Figure 4: Decoupling of the double occupied impurity-site of the bulk-chain at half-filling. 
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Figure 5: Impurity magnetization as a function of magnetic field for band filling 0.25 and several 
values of a. The dashed lines denote the asymptotic behaviour (4.20) near the critical magnetic field. 




Figure 6: Impurity magnetization as a function of magnetic field for band filling 0.9 and several 
values of a. The dashed lines denote the asymptotic behaviour for small magnetic field (see (4.15) 
and (4.17)) and near the critical magnetic field (4.20). 
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Figure 7: Impurity magnetization as a function of a for several constant magnetic fields and electron 
densities. 
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Figure 8: Number of electrons located at the impurity as a function of the bulk electron-density for 
fixed magnetic field = 0.1 and several values of a. The dotted line denotes the critical electron 
density (4.7) corresponding to H. 
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Figure 9: Number of electrons located at the impurity as a function of the bulk electron-density for 
various values of the magnetic field and a = 5. 
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Figure 10: The phase diagramm of the t-J model, shaded regions correspond to a ferromagnetic 
ground state. 
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Figure 11: Wilson-ratio as a function of the electron-density. 
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Figure 12: Derivative of the charge-stiffness with respect to the electron density at the crtitical 
electron density for a chain with 20 percent impurities as a function of a. 
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Figure 13: Charge-stiffness for a chain with 20 percent impurities and = 0.1 as a function of the 
electron density. The dashed hne denotes the stiffness at the critical electron density nc(a). 
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Figure 14: Charge-stiffness for a chain with 20 percent impurities and H = 0.5 as a function of the 
electron density. 




Figure 15: Charge-stiffness for a chain with 50 percent impurities and = 0.1 as a function of the 
electron density. 
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Figure 16: Spin-stifFness for a chain with 20 percent impurities as a function of the electron density. 



32 



